Emergence and control of heat current from strict zero thermal bias 
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It is an ever-growing challenge to develop nano-devices for controlling energy transport. An open question is 
whether we can create and control heat current at strict zero thermal bias, and if yes, how to do it. In this paper, 
we demonstrate that a nonlinear asymmetric system, when pushed out of equilibrium, can produce heat current 
in the absence of a thermal bias. The emergence and control of heat current over a broad range of parameters 
are studied. Our results reveal the following three necessary conditions: non-equilibrium source, symmetry 
breaking, and nonlinearity. We also demonstrate that when heat baths are correlated, symmetry breaking is 
sufficient to generate heat current. 
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Understanding heat transfer at the molecular level is of fun- 
damental and practical importance Oil. Recent years have wit- 
nessed a fast development in the emerging field of phononics 
121, wherein phonons, rather than an annoyance, can be used 
to carry and process information. To manipulate and control 
phonon transport (heat current) on the molecular level, various 
thermal devices 10, Ht] have been proposed. On the other hand, 
experimental works such as thermal rectifier j!4|| and nanotube 
phonon waveguide js]] have been carried out. These theoreti- 
cal and experimental works render the heat current to be con- 
trolled as flexibly as electric current in a foreseeable future. 

Heat transfers spontaneously from a high temperature to 
a low one; thus, the control of heat current has been so far 
based on the control of the temperature gradient. However, 
a large temperature gradient is essentially difficult to main- 
tain over small distance in practice, especially at nanoscale. 
Consequently, a natural question is raised: can we create and 
control heat current in the absence of (or against) thermal bias 
at nanoscale; if yes, then how do we do that? 

Inspired by ideas from Brownian motors |@], originally 
devised for particle transport, a few studies have revealed 
the possibilit y o f pumping heat against thermal gradients at 
nanoscale IIT HlOtl . A molecular model with modulated energy 
levels has been found to perform the heat pumping operation 
lUt]. While the microscopic oscillator system, though built on 
the similar principles, fails to perform the pumping l|9fl. Thus, 
it is still not clear what the requirements are for the system to 
show such functional effect. In this paper, we attempt to an- 
swer this novel and important question: how do we create and 
control heat current at strict zero thermal biasl 

It is noted that some interesting works reveal that nonzero 
heat cuiTent survives when one bath temperature is driven but 
with equal avera ge ( but different at any instant) to the other 
bath temperature lUOll . However, this reported behavior can be 
understood through the Landauer formula for the heat current 
d: J = / dioujT{Lo)[ri{Lo,TL) - vi^^Ta)], where T{uj) 
is the transmission coefficient and ri{uj,T]:^^]i) is the Bose- 
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FIG. 1: (Color online). Schematic illustration of one-dimensional 
two segment FK lattice, being coupled to two isothermal baths with 
oscillating temperatures = Tnit). 



Einstein distribution. Considering the temperatures /r are 
driven around the same average Tq, Taylor expansion gives: 
77(t^,TL)-77(c^,T«) ~ 77'('^,ro)[TL(t)-rfl(i)]+r/"[ri(t)- 
T7{(i)]^/2. After the periodic average, the first term vanishes 
while the second order survives which produces the nonzero 
current. 

Therefore, in a stark contract to the above proposals, we 
keep strict zero thermal bias at every instant through our stud- 
ies. This seems to be a small step, but it is a revolutionary one 
and has a completely different physics. It is under this strict 
zero thermal bias that our results uncover these three follow- 
ing conditions for the emergence of heat cuiTent at zero ther- 
mal bias: non-equilibrium source, symmetry breaking, and 
nonlinearity. Moreover, our simulation and analytic results 
reveal a phenomenon that symmetry breaking is already suffi- 
cient, if the two heat baths are correlated. 

Our system consists of two segment Frenkel-Kontorova 
(FK) chains llT2l ITsll coupled together by a harmonic spring 
with constant strength fci„t as depicted in Fig. [T] The Hamil- 
tonian can be written as: 



H = HL + ^{qN,,L-qi,Bf + Hn, (1) 

where the Hamiltonian of each FK segment reads: 
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FIG. 2: (Color online). Frequency resonance effect, (a) J versus 
Lu. (b) The cyclic local heat current Ji{t) (see definition in context) 
is illustrated at the optimal driving frequency a; = 27r x 10"'^. (c) 
Ti profile is shown for the same uj. N = 50, To = 0.09, AT = 
0.045, 20fc„t = 5/ci ^ kR = 1,5Vl = Fa = 5,7 = 0.5, 
throughout the paper, except specified. 



S stands for L or R, which represents the left or right segment 
with the same length. qi,s denotes the displacement from the 
equilibrium position for the ith atom in segment S and pi,s 
the corresponding momentum, a is the lattice constant, ks 
and Vs are the spring constant and the strength of the on-site 
potential of segment S. Two isothermal baths contacted with 
two ends are simulated by Langevin reservoirs with zero mean 
and variance (Ci/Af(OCi/Jv(^')) = 2^ kBTi^/RS{t — t'), where 
7 is the system-bath coupling strength. 

One question of interest is whether, due to the spatial asym- 
metry of the system, the thermal fluctuation in heat baths can 
induce a net heat current in a given direction. We argue that 
it is impossible. The situation would be a perpetual machine 
of the second kind extracting useful work out of ambient ther- 
mal reservoirs of vast energy surrounding us. Unfortunately, 
the second law of thermodynamics rules out the hiding place 
of the Maxwell demon, no matter how smart you design the 
system. It seems that any design to generate a heat current 
without thermal bias is foolish and even a quackery in the face 
of the second law. However, the second law works at thermal 
equilibrium only. 

In this paper, we drive the system out of equilibrium by pe- 
riodically oscillating two isothermal baths simultaneously, as 
TLit) = TR{t) = To + Arsgn(sinwt), where Tq is the ref- 
erence temperature. Under the time-varying heat baths, in the 
long-time limit, the local temperature of site i is time periodic, 
namely, T,{t) = mqf{t)/kB = T,{t + 2tt/u}). Similarly, flie 
time-dependent local heat current has the same periodicity: 
Ji{t) = kqi{t){q^{t) " qt+i{t)) = Ji{t + 2tt/lli). Therefore, 
within one period t G {0,2tt/uj), we can define the cyclic 
local heat current averaged over the ensemble of periods af- 
ter the transient time: Ji{t) = i X)fe=i "^^(^ + 2fc7r/cj) (n is 
the number of periods) as well as the cyclic local temperature 
^*(*) = Ti ELi + 2fc7r/cj). Thus, the net heat current 
and the effective local temperature read J ^ ^ J^'^ Ji{t)dt 
and Ti = ^ J^'^ Ti{t)dt, which are the same as the long- 
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FIG. 3: (color online). Temperature tuning effect (upper panel). 
UJ — 2iT X 10^^. The middle and bottom panels show Ti profiles and 
Ti{t) patterns with To = 0.07 (a) and 0.16 (b). At low temperature 
case (a), the left segment has lower conductivity which results in 
slower heat dissipating. Thus, more energies cumulate at the left 
part, making J from left to right. While at high temperature case (b), 
the right part has lower conductivity which induces the reversal of J. 



time average. 

For convenience of numerical calculations, we use dimen- 
sionless parameters by measuring positions in units of [a], 
momenta in units of [a{mkR)^/^], spring constants in units 
of [fcfl], frequencies in units of [(fc^/m)-'^/^], and tempera- 
tures in units of [a^fcij/fes]. A fixed boundary condition is 
applied and the equation of motion is integrated by the sym- 
plectic velocity Verlet algorithm with a time step of 0.005 for 
a sufficiently long time to guarantee the nonequilibrium stably 
periodic states. 

By adjusting the driving frequency lu of the two isothermal 
baths simultaneously, we obtain a nonzero heat current, which 
can be maximized at a moderate a;, as shown in Fig. [2^. It is 
intuitive that the heat current vanishes at both high and low 
frequency regimes. In the fast-oscillating limit lu ^ oo, ther- 
mal baths are driven so fast that two ends of the lattice cannot 
respond accordingly. The system only feels the same time- 
averaged temperature Tq at the ends, which yields J 0. 
In the adiabatic limit a; — !• 0, the system reduces to its equi- 
librium counterpart without thermal gradients so that J — 0. 
The emergence of heat current happens only when two seg- 
ments of the system respond differently. 

The Ji{t) pattern at the optimized frequency is shown in 
Fig. I2J3. At the first half driving period with positive tempera- 
ture variation, the heat transports from the two ends to the cen- 
tral part. While at the second half driving period with negative 
temperature variation, the direction of heat current reverses 
from the central part to the two ends. Eventually, the net heat 
cuiTent emerges due to the asymmetry of heat conduction re- 
sulting from the asymmetrical segment structure. Moreover, 
we illustrate the Xi profile in Fig. It indicates that more 
energies are cumulated at the left segment part which in turn 
induces the net heat current from left to right although two 
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FIG. 4: (Color online). System size effect. J versus cj with vari- 
ous sizes N for two different reference temperatures Tq: (a) 0.09; 
(b) 0.12, with AT = To/2, (c) Uc versus iV indicates the ballistic 
transport and the normal diffusion. 



heat baths hold the same temperature. As we mentioned, the 
emergence of heat current here does not break down the sec- 
ond law of thermodynamics since the system is pushed out 
of equilibrium. The only non-equilibrium source driving the 
heat current is the oscillating temperature TL{t) and Tfj{t) 
other than the thermal gradients. This phenomenon is some- 
how similar to that resonance in the particle current observed 
in the temperature ratchet ifTill . 

The reference temperature To is a very important param- 
eter since thermal conductivity is generally temperature de- 
pendent in many nonlinear lattices ifTsll . In the upper panel 
of Fig. |3] we show that the direction of the net heat current 
reverses as the reference temperature To increases and then 
saturates at high temperatures. To gain more insights into this 
reversal phenomenon, we depict the effective local tempera- 
ture profile Ti and the cyclic local temperature pattern Ti {t) 
with two typical values of Tq in the middle and lower panels 
of Fig. [3] It shows clearly that at low temperatures, energy 
dissipates faster at the right segment with high thermal con- 
ductivity which makes heat cumulated at the left part. Thus, 
there is a net heat current flowing from left to right. While 
at high temperatures, the scenario is reversed. In other words, 
the heat current flows from the segment with lower conductiv- 
ity to the higher one and the reversal of heat cuiTents results 
from the order reversal of thermal conductivities of two FK 
segments as temperature increases. Further, we check the am- 
plitude tuning effect by varying AT. As expected, it is found 
that the larger the AT, the larger the magnitude of J. 

By varying N, we find that the maximum value of J in- 
creases with system size and then saturates as shown in Fig. 
|4] Moreover, the optimum frequency coc decreases as N in- 
creases. This "redshift" can be understood from thermal re- 
sponse time ifioll . The heat conduction in the FK lattice fol- 
lows Fourier's law when the system size is larger enough lfl3ll . 
thus, the response time, r ^ N'^, characterizing the time scale 
for the energy to diffuse across the system. Therefore, the 




FIG. 5: (Color online), (a) Thermal bath correlation effect, c = 

2^kBTi^/ji{t). The upper black arrow marks J, without driving, un- 
der correlated baths with N = 8; while the bottom red arrow marks 
J in the adiabatic limit, with the same condition. Their discrepancy 
indicates the geometric phase induced heat current, (b) Correlation- 
induced J versus N. Results are calculated from the analytic Eq. 
^ and are checked that they are the same as the simulation results. 
AT = Vl = Vr = 0,c = 2jkBTo. 



characteristic frequency scales as cj^ c>c iV^^, which explains 
the redshift of the optimum frequency luc when N increases. 
More interestingly, we find that when N is decreased, the di- 
rection of J can be reversed and then the magnitude of J in- 
creases although negative. The optimum frequency ujc in this 
small size regime, scales as N^^ instead of the scaling N^'^ 
of the normal diffusion as shown in Fig. 2};. This is because at 
small N regime, the phonon transports ballistically. Since the 
right segment is more rigid than the left one (fc/j = bk^), the 
energy transports faster in the right part which induces J from 
right to left so as to reverse the direction. In other words, the 
crossover from normal diffusion to ballistic transport owing to 
the size reducing induces the reversal of the net heat current . 

In all simulations so far, two isothermal baths are indepen- 
dent, since (Ci(i)CAf (^')) = 0- Now we introduce the nonzero 
variance {£.i{t)^N{t')) = cS{t — t') to study the thermal bath 
coiTelation effect. We find that the coiTelation effect is sig- 
nificant only at small sizes (in which the energy transports 
ballistically) while fading away at large sizes in which energy 
transports diffusively, as shown in Fig. |5^. When = 8, 
we even obtain a nonzero heat current in the adiabatic limit 
a; — > 0. It implies that there may exist a"geometric phase" 
in the dissipative and stochastic system |16;] we used here. 
The geometric phase results from the nonzero area enclosed 
by periodic variation of parameters in parameter space. Even 
when the driving is extremely slow, approximating to equilib- 
rium state at every instant, this geometric phase-induced heat 
current still survives. This is an interesting topic deserving 
further investigation. 

Moreover, in contrast to the nonlinear FK lattice, we find 
surprisingly that, for a pure harmonic system, even in the ab- 
sence of external driving, and without thermal bias, the heat 
current emerges when the two thermal baths are correlated. To 
understand this correlation-induced heat current, we formal- 
ize the heat conduction of harmonic systems by the Rieder- 
Lebowitz-Lieb method UtIi and have the steady-state equa- 
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tion: 

A-B + B- ^b, (3) 

where A = \ - , Z) = „ ^ ■ -^ h" = 

kim — kij is the Laplacian matrix where kij de- 
notes the spring constant between adjacent oscillator i and 
j. Tin ~ li5im{n = 1,-/V) is the dissipation matrix. 
Cij = 2j.ikBTi{6iiS ji + SinSjn) + ciSiiSjN + SinSji) and 
the second term depicts the correlation effect. The solution 

B of Eq. ^ has four blocks: B = ( f f f ) , where 

V Bqp Bpp J 

{Bqp)ij ~ {qiPj) denotes the position-velocity correlation re- 
lating to J. For two oscillator system, it is easy to obtain the 
expUcit solution of heat currents: 

J _ kint[{kL - kR)c + 2-f{TL - Tji)kint] , . , 

" {kL - fc/?)2 + 272 (fcL + kji) + 4^^k„,t + 4fcf„, ■ ^ ' 

It shows clearly that even without thermal bias (T^ = Tr), 
the correlation c still can induce nonzero J in the asymmetry 
harmonic chain (fc^ 7^ ku). This nontrivial correlation effect, 
represented by the nonzero off-diagonal of Cij, is reminis- 
cent of the quantum coherence ||3ll9!], which provides the 
off-equilibrium source to create heat currents. More interest- 
ingly, we find that when N increases the correlation-induced 
heat current oscillates and reverses its direction periodically, 
as shown in Fig. [Sji- 

In summary, we have demonstrated the emergence of heat 
currents and their direction control in two segments FK lattice 
without thermal bias. By periodically oscillating two isother- 
mal baths simultaneously, we can obtain a maximum heat cur- 
rent at an optimum driving frequency. The resonance effect 
with respect to other parameters has been studied systemati- 
cally as well. We have found that the direction of the emerg- 
ing heat current can be reversed by either tuning the reference 



temperature or tuning the system size, and the optimum fre- 
quency can be decreased by increasing the system size. Our 
results reveal the following three necessary conditions for the 
emergence of heat current without thermal bias: (1) nonequi- 
librium source, which is induced by the periodically oscillat- 
ing baths although isothermal, and in turn breaks the underly- 
ing detailed balance. Also, we can periodically drive other pa- 
rameters such as kint to generate the nonequilibrium source so 
as to create heat current (not shown); (2) symmetry breaking, 
which results from the two asymmetric segments construction 
defining a preferential directionality; (3) nonlinearity, which 
comes from the on-site sinusoidal potential in the FK model. 
In fact, when V[^/fj = 0, the system reduces to an asymmet- 
ric harmonic chain and we find that the heat current vanishes 
(not shown). However, if the correlation between two baths 
is introduced, symmetry breaking is sufficient to create heat 
current. 

We should point out that the model proposed here might 
be realized experimentally. For a typical atom, a ~ 1 
A, m ^ 10^2^ — 10^2^ kg, which yields the frequency 
unit [{kn/my/^] ~ 10^'^ s^^ and the temperature unit 
[a^kR/ks] ^ 10^ - 10"* K. The typical value of Tq = 0.1 and 
u! 10^'^ — 10~^ corresponds to the physical temperature 
Tr - 10^ - 10^ K and the physical frequency uj^ ^ 10^ - 10^ 
MHz, which is in the ultrasonic and microwave regimes. 
Thanks to the redshift effect, we are able to obtain lower opti- 
mum driving frequencies in practice by enlarging the system 
size. The thermal bath correlation might be implemented by 
imposing some common external thermal noise, or introduc- 
ing entanglements between two heat baths lfl9ll . We hope the 
present study will stimulate experimentalists to search possi- 
ble realizations and technological utilizations. 
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